Pauli-limited upper critical field in dirty <i-wave superconductors 
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We calculate the Pauli-limited upper critical field and the Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) instability for dirty d-wave superconductors within the quasiclassical theory using the 
self-consistent t-matrix approximation for impurities. We find that the phase diagram depends 
sensitively on the scattering rate and phase shift of nonmagnetic impurities. The transition into 
the superconducting state is always second order for weak (Born) scattering, while in the unitarity 
(strong) scattering limit a first-order transition into both uniform and spatially modulated super- 
conducting states is stabilized. Contrary to general belief, we find that the FFLO phase is robust 
against disorder and survives impurity scattering equivalent to a T c suppression of roughly 40%. Our 
results bear on the search of FFLO states in heavy-fermion and layered organic superconductors. 

PACS numbers: 74.25.Ha, 74.25.Dw, 74.81.-g, 74.25. Op 
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Introduction. In typc-II singlet superconductors a 
magnetic field suppresses superconductivity for two rea- 
sons: (1) the phase of the Cooper pair wave function 
couples to the vector potential resulting in the appear- 
ance of vortices; (2) Zeeman coupling of the magnetic 
field to the electron spins polarizes and splits the con- 
duction band, which destroys superconductivity when 
the loss in magnetic energy equals the energy gain from 
pair condensation^ 0, [H, 3 ■ This latter mechanism 
is referred to as Pauli limiting and leads to a first or 
second-order transition from the normal (N) to supercon- 
ducting (SC) state depending on the value of the mag- 
netic field. It has been predicted that a clean system at 
high fields can remain superconducting beyond the Pauli 
limit by forming the nonuniform Fuldc-Fcrrcll-Larkin- 
Ovchinnikov (FFLO) state with a spatially modulated 
order parameter This state, however, is suppressed 
by disorder Q. 

In contrast to conventional (isotropic s-wave) su- 
perconductors, unconventional (d-wave) superconductors 
are affected by nonmagnetic impurities even at zero field; 
scattering averages the gap over the Fermi surface and 
suppresses T c . The different rates of suppression of the 
uniform and FFLO states determine the phase diagram 
in the field-temperature {B-T) plane. Agterberg and 
Yang found that in two-dimensional (2D) <i-wave su- 
perconductors with purely Zeeman coupling, the N-SC 
transition is of second order at all T, with the Larkin- 
Ovchinnikov (LO) modulation, A^o ~ cosq-R, and the 
uniform (USC) state, Ausc = const, favored at low and 
high temperatures, respectively, and a narrow intermedi- 
ate T region, where the nodclcss Fuldc-Fcrrcll (FF) state, 
App ~ e 4q R , is stabilized. Ref. || reported that under 
combined orbital and Zeeman coupling in impure d-wave 
superconductors the first-order transition into the vor- 
tex state appears at intermediate temperatures. Very re- 
cently, Houzet and Mineev 0] studied orbital and impu- 
rity effects in s-wave and d-wave Pauli-limited supercon- 



ductors and concluded that orbital effects are necessary 
for a first-order transition to occur in 2D d-wave super- 
conductors. In contrast, for s-wave in 3D the transition 
to the FFLO state is first-order [lol[ll|. 

Remarkably, our understanding of impurity effects in 
nonuniform states is still incomplete. In Refs. 0, @, @| 
the discussion was limited to weak (Born) impurity scat- 
tering and focused only on the Ginzburg-Landau (GL) 
regime close to the onset of the FFLO instability, using 
an expansion in the modulation wave vector q. However, 
q = |q| increases rapidly to values comparable to the in- 
verse superconducting coherence length, q£o ~ 1, so this 
expansion quickly becomes invalid away from the critical 
point. 

In this Letter, we present a microscopic treatment of 
impurity effects on the superconducting states in purely 
Pauli-limited quasi-2D d-wave superconductors. Impu- 
rities are treated in the self-consistent i-matrix approx- 
imation (SCTA) covering the weak (Born) and strong 
(unitarity) scattering limits [t|. The latter limit, never 
considered previously, is especially important because of 
a search for FFLO-like states in heavy-fermion and lay- 
ered organic superconductors [l3j], where impurity scat- 
tering is strong Our approach is not limited to an 
expansion in q, and hence is valid for any temperature 
and impurity concentration along the second-order upper 
critical field B C 2- We show that the phase diagram of a 
Pauli-limited dirty d-wave superconductor is very differ- 
ent for nonmagnetic impurities in the Born and unitarity 
limits. The differences originate from the dependence 
on scattering strength of quartic and higher order coeffi- 
cients in the GL functional. The first order N-SC tran- 
sition, absent for Born scattering, is stabilized by strong 
impurities, and is therefore expected in heavy fermion 
systems. 

Quasiclassical equations. We follow Refs. l3 : EH and 
solve the quasiclassical equations for the 4 x 4-matrix 
Green's functions in particle-hole and spin space, which 
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satisfy the normalization condition, g 2 = —tt 2 1, and the 
transport equation, 

[ie m r 3 - nB ■ S - A(R, p) - a" np (R; e m ) , ( 1 ) 

g(R, p;e m )] + ihv f (p) ■ V R 5(R,p;e m ) = 0. 

Here /z is the magnetic moment, e m — TtkBT(2n + 1) 
are the Matsubara frequencies, A is the mean-field su- 
perconducting order parameter depending on the coor- 
dinate, R, and momentum direction, p, at the Fermi 
surface with velocity vj. The electron spin operator is 
S = er^(l + T3) + cr*2(l — T3). The Pauli matrices er and 
r operate in spin and particle-hole space, respectively. 
Eq. |T]) is complemented by self-consistency equations for 
A and the impurity self-energy a lmp . We use h = = 1- 

In the SCTA a lmp = ni mp t, with impurity concentra- 
tion riimp. For isotropic scattering the t- matrix satisfies 
t(R; e m ) = u Q l + u A/"/(#(R, p; e))p t(R; e m ), where an- 
gular brackets (...) denote a normalized Fermi surface 
average. The strength of the nonmagnetic impurity po- 
tential, uq, is expressed via the isotropic scattering phase 
shift, = arctan(7ruoA//); A/ is the density of states per 
spin at the Fermi surface. For Born (unitarity) scattering 
8q = (Sq = 7r/2) and the normal-state scattering rate 
T = l/2rAr = T u sin 2 S , with T u = n imp /iiNf. 

If we choose the direction of the spin quantization 
along B = Bz (which is allowed if the hamiltonian has 
spin-rotation symmetry in the absence of the field), both 
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g and a lmp have block-diagonal structure corresponding 
to the two spin projections. Hence, the quasiclassical 
equations for the spin-up and spin-down sectors decou- 
ple , and we solve separately for the diagonal, g s , and 
off-diagonal, f s , f' s , components of <?, with s = ±1{|, I}, 
with the constraint g 2 — / s / s ' = — tt 2 . However, both spin 
projections enter the self-consistency equation for A. We 
assume a separable pairing interaction y(p)y(p'), where 
3^(p) gives the angular dependence of the gap function 
with the normalization (y 2 (p)} = 1. For A(R, p) = 
A(R)y(p), we find 

A(R)ln|-=T^^(R,p; £m ))--^^) , (2) 

-imp _ o ( COt5 Q + (g s )/TT (/ S )/7T \ / o\ 

5 " *V a') A cotS - (gj/n ) ' [6) 

Here J"(R,p;e TO ) = ^(p)[/ T (R, p; e m ) + / X (R, p; e m )\ 
and S s = r/[l - Tr^sin 2 S ((g a ) 2 - (/.)</.') + tt 2 )]. To 
calculate the B-T phase diagram, we derive the 
Ginzburg-Landau functional (expansion in A for arbi- 
trary q) by taking A(R) = A q exp(iq-R) and solving 
Eqs. dTJ- ((3j) together with the normalization condition 
for g to third order in A. We substitute the n-th or- 
der solutions / s (1) , / s (3) into Eq. © to obtain the GL free 
energy difference between the SC and N states, 



An GL = ^a(T, 5; q)|A q | 2 + £ -f3(T, B; qi, q 2 ; q 3 , q 4 ) A qi A q2 A q3 A q4 <5 qi+q2 , q3+q4 , (4a) 
q qiq2q3q4 

a(T, B;q) = In - 2irT £ Re ((^D" : ) - e" 1 ) , (4b) 

/3(T, J B;q 1 ,q 2 ;q 3) q4) = ttT £ Re ( /^q^^ ^+^+W*>/ 4 \ - rY qiq2q3q4 1 , (4c) 

t - (- -sin 2 <5n^ fe (2) e (2) +e (2) e (2) "l -sin 2 <5 n fe (1) e (1) e (2) 

1 qiq2q 3 q4— l 2 bm °0 I ^qiqa^qiqi + U qiq4 U q2qs j hlU °° ^qi ^q 3 U q 2 q4 

+0(1)0(1)0(2) +0(1)0(1)0(2) +0(1)0(1)0(2) _ 2 e (1) e (1) e (1) e (1)N l (4d) 

+ ^qi ^q 4 U q 2 q 3 + ^q 2 U q 3 ^qiq4 + ^q 2 ^q4 U qiq 3 ZU qi U q 2 ^q 3 ^q 4 J ' (* U ) 



where we introduced the angular averages O qi = 

(y^D^, Oqfq. = (j^q^q^q/), and defined 

r /q = h w f ' 1' and = £m + r + ^O"- 8 + r ?q)- Wc 
introduced j) q = ^ + ^, q with X,q = re q 1} . 

Results. The second-order N-SC transition is deter- 
mined from the GL coefficient a(T, B;q) =0 and de- 
pends only on T, but not on i5o- Thus the transition 
line is independent of the phase shift. An instability 



into the modulated FFLO state becomes possible below 
T q ^o, where the maximal B c i is found for q =/= 0. This 
occurs when the GL coefficient k in the g-expansion of 
a(T, B; q) ps ao + nq 2 becomes negative, 

For d-wave SC the last term vanishes, since (3^?q) = 0. 
In contrast, the quartic term in the GL functional ex- 
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FIG. 1: (Color online) The critical fields of 
Pauli-limited dirty d-wave superconductors 
for transitions to USC (first-order, dashed) 
and FFLO states (second-order, solid) for 
q along nodes (squares) and antinodes (cir- 
cles) . Superconductivity sets in at the high- 
est B C 2- (a) Pure case: Below Tp — T q ^o 
the N-USC transition is below the second- 
order N-FFLO transition, (b-d) Born and 
unitarity impurities split Tp and T q ^o , thus 
Tp < Tq^o < Tp. A modulation along 
antinodes quickly suppresses B C 2 relative to 
q||node, panel (b). For Born impurities the 
second-order FFLO transition is above the 
N-USC line at all temperatures, as for the 
clean case. In the unitarity limit the first- 
order N-USC transition preempts a N-LO 
transition, see panel (d). 
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plicitly depends on the scattering phase shift, Eqs. (j4c|l - 
(|4dj) . For example, it controls the location of the first- 
order transition to the USC state, Tp, which competes 
with the FFLO instability. In unconventional supercon- 
ductors (y) = 3-1,0 = and the critical point Tp is deter- 
mined by a sign change of the GL coefficient (3 at q = 0, 

&=^i;m- T -^r^)- TO 

For r = 0, both k and (3q become negative at exactly the 
same temperature, T q ^o = Tp ~ 0.5615T c o. Since the 
transition into the FFLO state has a higher critical field 
at any temperature T < Tp. the first-order transition is 
superseded by the onset of the FFLO state [l(| ■ 

A comparison of k and (3q shows that in dirty uncon- 
ventional superconductors T q ^o = Tp only for Sq = n/4. 
For Born (B) and unitarity (u) scattering (3o depends 
on Sq, such that Tp and Tp 1 shift in opposite directions 
relative to Tq^ , hence Tp < T q ^o < Tp as shown in 
Figffl The latter inequality is especially important since 
it shows that for strong scatterers Pauli limiting leads to a 
first-order transition into the USC state at high fields/low 
temperatures in the B-T phase diagram. As the sys- 
tem becomes dirtier, i.e., the lifetime decreases, these 
characteristic temperatures are suppressed to zero in the 
following order, T® -> at T/irT cQ > 0.18, T q#0 -> 
at T/ttTcq > 0.20, and T£ -> for T/ttT c0 > 0.22. Note 
that for larger T the N-USC transition line is of second 
order at all T. 

Fig. Q] gives the upper critical field lines for different 
states. Second-order transition lines are found by the 
largest spatial modulation vector q = Q that maximize 
B c2 . In clean d-wavc SC [IE [13, 0, El the modulation 
is along a gap maximum (antinode) at low T/T c q < 0.06, 
and along a gap node for 0.06 < T/T c0 < 0.56, see Fig. 
[TJa). However, already for small impurity scattering, 



T/ttT c o > 0.02, the critical field for q\\antinode is low- 
ered below B^ node j anc j t ne stable configuration is with 
q||node over the entire range of existence of the FFLO 
state, see Fig. [IJb). 

Determining the first-order transition lines of B c2 re- 
quires a self-consistent calculation of the full free en- 
ergy functional, the details of which will be given else- 
where ■ We find that in the Born limit the first-order 
transition is always below B^ 2 FLO , in agreement with 
0,111 . In contrast, in the unitarity limit T q ^o < Tp and 
B^ 2 FLO is below the first-order transition to the USC 
state, see Figs. QJb-d). 

For intermediate impurity scattering, the phase dia- 
gram is given in Fig. [2j To determine the structure 
of the SC state near B c2 , we analyze the GL free en- 
ergy, Eq.(j4]), for four possible phases: USC [A(R) = 
Aj/sc], FF with a single Fourier component Qi = 
(Q,0) [A(R) = A FF exp(iQx)}, LO with {Qi,Q 3 } = 
{(±Q,0)} [A(R) = Alo2cosQx], and square lat- 
tice (SQ) with {Qi,Q 3 ,Q2,Q 4 } = {(±Q,0),(0,±Q)} 
[A(R) = Asq2(cos Qx + cos Qy)]. The x,y— axes are 
along the gap nodes. For each phase, we calculate 
Anf L = -a 2 //3 z , with Pff = Willi, Plo = Pun + 
2/? 1313 and (3 SQ = 0.5(Aui + W1212 + Wiais + Wuu + 
2/3i324), where /3 ijk i = P(T, B; Q 4 , Qj-; Q fe , Qz). Along 
the second-order transition line the phase with the low- 
est positive value of (3 has the lowest energy. 

For Born impurities (Fig. [2] right), j3i > for all 
nonuniform states, and the LO state is favored in most 
of the phase diagram except a small region below Xq^o, 
where the FF phase is stabilized for the impure case 
0, Q . Analysis of Afl GL indicates that this phase is sep- 
arated by a second-order transition from the USC and by 
first-order from the LO state. 

The situation is very different for strong impurities 
(Fig. |2 left). Following the B c2 line from T C (B = 0) to 
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FIG. 2: (Color online) The phase diagram for T /tvT c0 = 0.14. 
Left panel: transition into the LO state at low T becomes 
first order above T[q I! for unitarity impurities. There is also 
a region of a first-order transition into the USC state below 
Tp. Right panel: Born impurities result in second-order tran- 
sitions. The LO state is favored in large parts of the phase 
diagram over the FF state, except near T q7 £o 0, H ■ The tran- 
sition line shown between the LO and USC states (dotted 
line) is qualitative. 
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FIG. 3: (Color online) Phase transitions along the B C 2 line in 
the T-T c plane. The (black) solid line is T C (T) at B = 0. For 
Born impurities the LO state exists for T < Tq^o, with a small 
region occupied by the FF state. For unitarity scattering a 
first-order transition appears into the USC and LO states. 
At low T a square lattice FFLO state (SQ) [3 is rapidly 
suppressed with increasing F [22l |. 



lower T, we reach the critical point Tp, below which the 
N-USC transition is of first order. At T — > the transi- 
tion is second order into the LO state, but becomes first 



order above T 
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We estimate where the first-order 



N-USC and N-LO lines meet. However, determining the 
location of this point and the LO-USC transition line re- 
quires a fully self-consistent treatment of the nonuniform 
problem 15j, [21| , which is beyond the scope of this work. 

Conclusions We summarize our results in Fig. [3l 
where we show all states that arise along the upper crit- 
ical field line for fixed T. For nonuniform states, we 
only consider modulations along gap nodes, since states 
with q\\antinode arc destabilized even faster by impuri- 
ties. We find for dirty e?-wave superconductors that the 



FFLO state is quite robust and survives impurity scatter- 
ing equivalent to ~ 40% of T c suppression or a mean-free 
path £ of £o/^ ^ 0.16. This result is important for the 
search of an FFLO state in doped Cc-115 [HI], an d other 
heavy-fermion and layered organic superconductors. 

Notably, the differences between weak and strong im- 
purity scattering are significant. In the Born limit Tp 
is suppressed below the onset of the nonuniform state, 
Tp < Tq^o, and the transition is always of second order. 
Impurities stabilize a narrow region of the Fulde-Ferrell 
state just below Tq^o- In contrast in the unitarity limit 
(relevant to recent experiments) T q ^o < Tp and the first- 
order transition into the uniform state preempts a mod- 
ulated state. Below T ~ T q ^ the transition into the 
Larkin-Ovchinnikov state begins as a first-order line and 
becomes second-order at lower T. Importantly, in this 
limit the interplay of Zeeman splitting and disorder, even 
without orbital effects, drives the transition between the 
normal and superconducting state first order. 
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